Nonlinear Meissner Effect in CuO Superconductors 
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Recent theories of the NMR in the CuO superconductors are based on a spin-singlet d x 2_ y 2 order 
parameter. Since this state has nodal lines on the Fermi surface, nonlinear effects associated with 
low-energy quasiparticles become important, particularly at low temperatures. We show that the 
field- dependence of the supercurrent, below the nucleation field for vortices, can be used to locate 
the positions of the nodal lines of an unconventional gap in momentum space, and hence test the 
proposed d x 2_ y 2 state. 
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C***'* ■ Much of the interest in the CuO high T c superconductors is due to the belief by many theoreticians that these 
materials represent new, or at least novel, forms of superconductivity. Theories based on very different excitations 
{e.g. spinless fermions, anyons, and traditional charged fcrmions) have been proposed to explain the normal and 
superconducting properties of these materials. Superconductivity from anyons has received considerable attention. 
I . Similarly, several theories based on BCS-like pairing, either of novel excitations or of novel pair states of conventional 
' charged fermions have been proposed [ Q]. All of the theories of the superconducting state, except the conventional 
BCS pairing state, have the property that they break at least one additional symmetry of the normal state besides 
gauge symmetry. While there has been little compelling evidence to identify the order parameter of the high T c 
superconductors with any of these unconventional states, recent theories, based largely on interpretations of NMR in 
both the normal and superconducting states of YBCO, have focussed on the spin-singlet, d x 2_ y 2 state as the order 
— t \ parameter for these CuO superconductors [ J^]]. In this letter we show that the field dependence of the supercurrent 
can be used to locate the positions of the nodal lines of the gap in momentum space; we suggest new experiments to 
' resolve the gap structure of the d x 2_ y 2 state, or other unconventional BCS states. 
""j" 1 , Based on ideas from Hammel, et al. and Shastry [ 0]], a semi-phenomenological model of a nearly AFM Fermi liquid 
that accounts for the anomalous NMR properties of the CuO superconductors has been proposed [ am]. Monthoux, 
et al. [ O] recently argued that the same AFM spin fluctuations that are responsible for the anomalous normal state 
magnetic properties are also responsible for superconductivity in the layered CuO systems. These authors examine 
the weak-coupling gap equation with a pairing interaction given by the exchange of a single spin-fluctuation, and find 
that the superconducting state with d x 2_ y 2 symmetry is the preferred solution. 

Several authors [ ||] have extended the nearly AFM Fermi liquid into the superconducting state, and have concluded 
that the Knight shift and anisotropy of the spin- lattice relaxation rates in YBCO are most reasonably fit with a singlet, 
d x 2_ y 2 order parameter. Since the arguments put forward in favor of a d x 2_ y 2 order parameter are connected with 
the spin excitation spectrum, independent tests of the d x 2_ y 2 state based on the orbital structure of the gap or its 
excitation spectrum are needed to identify the order parameter. 

The important feature of the quasiparticle excitation spectrum, which is required by symmetry for a clean super- 
conductor with a d x 2_ y 2 order parameter, is that the excitation gap vanishes along four lines ('nodal lines') located 
at the postions | fc^ | = \k y \ in the 2D plane and running along the length of the Fermi tube [see cross-section in Fig. 
1]. These nodal lines imply low-energy excitations, at all temperatures, that give rise to power law deviations of the 
penetration depth from its zero-temperature limit; i.e. 1/A(0) — 1/A(T) ~ (T/T c ) for T <C T c and a gap with line 
nodes, while deviations proportional to higher powers are obtained from gaps with point nodes [ 0]. At present 
there is no consensus 'for' or 'against' a true gap in the spectrum from penetration depth measurements. Several 
measurements of the penetration depth in CuO superconductors are interpreted as supporting the predictions of BCS 
theory with a true gap, while the analysis by Annett, et al. [ [0] of those same measurements [see also [ |8)] for original 
references] suggests deviations from a true gap, but gives no indication of a linear temperature dependence that would 
be consistent with a d x 2_ y 2 state. 

The main idea of this paper is that the field-dependence of the supercurrent may be used to locate the positions of 
the nodal lines (or points) of an unconventional gap in momentum space. This is a stronger test of the symmetry of the 
superconducting state than power law behavior for the penetration length. This gap spectroscopy is possible at low 
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temperatures, T <C T c , and is based on features which are intrinsic to nearly all possible unconventional BCS states 
in tetragonal or orthorhombic structures. The constitutive equation relating the supercurrent, j s , to the velocity, 
v s = \{d<j) + —A) (this definition avoids the introduction of an arbitrary mass scale), is linear only for velocities 
small compared with the critical velocity, v c — A/vf. Since the unconventional gaps of interest have nodal lines on 
the Fermi surface [ Q], nonlinear corrections to the constitutive equations become important, particularly at low 
temperatures. The low-energy states near the nodes of the gap lead to a quasiparticle contribution to the current that 
(i) persists at zero temperature, (ii) is a nonanalytic function of the velocity and (iii) exhibits a large anisotropy at 
low temperature that is determined by the positions of the nodes on the Fermi surface. We discuss the theory of the 
nonlinear current/ velocity relation at low temperature, and then show how the nonlinearity in the current leads to a 
novel basal plane anisotropy of the screening length and magnetization energy (or magnetic torque), for the d x 2_ y 2 
state when the field is parallel to the CuO planes. 

Consider a superconductor-vacuum interface which is parallel to the magnetic field, and assume the penetration 
length is large compared to the in-plane coherence length, A £. The important spatial variations are then contained 
in the local value for the superfluid velocity, which is effectively uniform on the scale of £. The order parameter 
achieves its equilibrium value for the local value of the velocity, and the resulting current is [ p0|], 

I = -eN f J d 2 s n(s) v f (s) ^a v (s) + jf d? [f(E + a v (s)) - f(E - <7„ (*))]} (1) 
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where a v — Vf(s) ■ v s is the shift in the quasiparticle energy due to the superflow, Nf is the total density of states at 
the Fermi level, n(s) is the angle-resolved density of states normalized to unity, Vf(s) is the quasiparticle velocity at 
the point s on the Fermi surface, and E = ^/^ 2 + |A(s)| 2 . The first line of eq.(Q) is useful for calculating the current 
in the low-temperature limit. The first term represents the supercurrent obtained from an unperturbed condensate, 
i.e. j s — —epv s , where p is the total density (p — \NfV 2 for a cylindrical Fermi surface). The term involving the 
Fermi functions is the current from the equilibrium population of quasiparticles. This formula shows clearly that the 
quasiparticles do not contribute to the current at T — for a superconductor with a true gap over the entire Fermi 
surface, provided the velocity is below the critical velocity, A m i n /vf. At finite temperature, the excitation term is 
responsible for the reduction in the superfluid density from p to p s (T) in the linear response limit. In addition, the 
quasiparticles give rise to a nonlinear reduction in the supercurrent. 

In order to motivate the proposal for detecting the nodal structure of the gap in an unconventional superconductor 
consider the effect of a magnetic field on the penetration depth in a conventional, isotropic superconductor. For a 
field parallel to the surface, the screening current is proportional to the applied field, i.e. j s (0) ~ cH/X. This linear 
relation breaks down because of pairbreaking induced by the velocity field, 



j a = p s {T) v s [1 - a(T) 
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for (v s /v c ) <C 1, where a(T) > and v c — A(T)/u/ is the critical velocity. One important point to note is that 
a(T) — > for T — ► 0, since the quasiparticle occupation is zero for all quasiparticles states with velocity below 
the critical velocity. The importance of the nonlinear current /velocity constitutive equation to the penetration of 
magnetic fields into a superconductor is qualitatively clear. The reduction of the current by the pair-breaking effect 
reduces the effective superfluid density, and, therefore, increases the effective penetration length. Since the current is 
proportional to the field in linear order, the correction to the effective penetration length is quadratic in the surface 
field [ 111)], 
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where H (T) = (3/4)(c v c {T)/e X(T)) is of order the thermodynamic critical field. 

Now consider the d x i_ y i order parameter, A(s) = A D (fc^ — k 2 ). The nodes in the gap imply that there is a 
quasiparticle contribution to the supercurrent even at T = 0, which is easily calculated in terms of the phase space of 
occupied quasiparticle states. Consider the case where the velocity is directed along the nodal line k x = k y \ v s = v s x' 
as shown in Fig. 1. For any non-zero v s there is a wedge of occupied states near the node opposite to the flow velocity. 
For a cylindrical Fermi surface the quasiparticle current is 
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/Ct c 
— cos 2 (a) y/ sin 2 (a c ) - sin 2 (a) , (4) 
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where a is the angle measured relative to the node at —x' and a c — sirT 1 (-s^ 2 -) is the maximum angle for which the 
quasiparticle states near the node are occupied. To leading order in (-r^- 1 ) we obtain a total supercurrent of 

7._ ( _ii^) M l_J*L- } , (5 ) 

for v s directed along any of the four nodes. Note that the current is parallel to the velocity, and that the occupied 
quasiparticle states reduce the supercurrent, as expected. Also the nonlinear correction is quadratic rather than cubic, 
as is obtained for the conventional gap, and with the characteristic scale determined by v = 2A /v/. Furthermore, 
the quasiparticle contribution to the current is a nonanalytic function of v s [ ]T^]]. 

Unlike the linear response current, the nonlinear quasiparticle current is anisotropic in the basal plane. A velocity 
field directed along the maximum direction of the gap, v s — v s x, produces two groups of occupied states, albeit with 
reduced populations because the projection of v s along the nodal lines is reduced by l/v2- The resulting current is 
easily calculated to be 

?. = (-|iV,.?)ft{l-^^L}, (6) 

which is again parallel to the velocity and has a quadratic nonlinear correction. However, the magnitude of the 
nonlinear term is reduced by l/\2. This anisotropy is due to the relative positions of the nodal lines and is insensitive 
to the anisotropy of the the Fermi surface or Fermi velocity because the quasiparticle states that contribute to the 
current, for either orientation of the velocity, are located in a narrow angle, a < a c ~ (v s /v ) <C 1, near the nodal 
lines. Thus, the occupied quasiparticle states near the nodes have essentially the same Fermi velocity and density of 
states; only the relative occupation of the states is modified by changing the direction of the velocity. 

The dependence of the supercurrent on the positions of the nodal lines in momentum space suggests that the 
anisotropy can be used to distinguish different unconventional gaps with nodes located in different directions in 
momentum space. For example, the d xy state, A ~ k x k y , would also exhibit a four-fold anisotropy, but the nodal 
lines are rotated by 7r/4 relative to those of the d x i_ y i state, while the state A ~ k x k y (k 2 — k y ) corresponding to the 
-A-2 g representation would exhibit an 8— fold anisotropy. 

This anisotropy in the current implies a similar anisotropy in the field dependence of the penetration length, which 
can be calculated from eqs.(||)-(||) and Maxwell's equation; in the gauge, d ■ v s = 0, and for surface fields parallel to 
a nodal line the equations reduce to, 

d 2 v s 47re 2 v s f \v s \\ 

Js[Vs\ = i 1 ~ — ^ • (7) 
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The magnetic field in the superconductor, b = (c/e)d x v s , is also parallel to a node, and has magnitude b = 
(c/e)(dv s /dz) with 6(0) = H. For field penetration into a bulk superconductor we define the effective screening 
length in terms of the surface impedance, 1/A e // = —(l/H)(db/dz)\ z= Q. We obtain 

1 1 / 2 H \ 1 1 / 1 2 H \ -„ 

T 1 - oTT - #11 node ; 1 = T 1 - -7k oTT ' H W antmode , (8) 



A e // A \ 3 H J ' A e // A \ -y/2 3 H a 

to leading order in H/H , where H = (v / X)(c/e) ~ </> /(A£) ~ H c . Observation of this anisotropy and linear field 
dependence would be strong indication of a d x 2_ y 2 order parameter. 

Sridhar et al. [ [O] have measured the field dependence of the penetration length for H parallel to the a and b axes of 
single crystals of YBCO. These authors report an in-plane penetration depth which obeys X(H, T) = A(T) + k(T)H 2 . 
At intermediate temperatures, T ~ 80 K, they measure k(80) ~ 1 A/G 2 , which drops by three orders of magnitude by 
T ~ 10 K; k(4) ~ 10~ 3 A/G 2 . The simplest interpretation is that the data is in good agreement with a full gap on the 
Fermi surface, in which case the d x 2_ y i state is not the order parameter of YBCO. However, we can ask if the data, 
so far, forces this conclusion. First of all, the d x 2_ y 2 state will also show a quadratic field correction to the effective 
penetration depth at intermediate and high temperatures, with a coefficient that is similar in magnitude to that of 
the conventional BCS prediction, and also drops rapidly with temperature. The question is, 'at what temperatures 
and fields does the effective penetration depth exhibit a linear field dependence?' The linear field regime is inferred 
from eq. ([[]) . The linear behavior arises from the nonanalytic dependence of j s on v s at T = for the nodal directions 
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where A(s) = 0. Non-zero temperature removes the singularity. For sufficiently small fields an expansion in v s 
is valid; however, the cross-over field is determined when the smallest Matsubara frequency becomes comparable 
to the superflow kinetic energy shift per particle, i.e. ttT ~ vtv s ~ 2A (H/H ). This implies a cross-over field, 
H x ~ H (T/T C ). Thus, to obtain a significant linear field regime down to H ~ H C \/1Q ~ 25 G requires temperatures 
T < 2 K. Thus, the experiments of ref.[ [^3[], which are reported down to 10 if, are apparently at the high end of 
the temperature regime where a linear term might become observable. 

We also consider the effect of impurity scattering on the linear field dependence of the effective penetration length. 
Quasiparticle scattering from impurities plays a similar role to that of temperature; the impurity lifetime removes the 
singular behavior of the current at sufficiently low velocities. The cross-over field due to impurity scattering can be 
obtained by including the impurity scattering self energy. The resulting current is given by eq. (Q) with the replacement 
of + (e n /2r) < l/y/il + |A(s)| 2 >, where r is the s-wave scattering rate in the Born approximation, 

and the average is over the Fermi surface. At zero temperature, impurity scattering leads to a new low energy scale [ 
JL4|l ]) e n = 4A e~ 7rrA °, which smooths out the linear behavior of the effective penetration length below the crossover 
field, H T = 2i/ e -7rrA °. This condition is considerably less restrictive than that for the temperature. For YBCO 
with a mean-free path of Z/£ = ittA ~ 10 the crossover field is of order 2H a e~ w <C H c ±. Thus, there should be no 
difficulty from impurity scattering in observing the linear behavior in the effective penetration length for moderately 
clean YBCO crystals. 

Another test of the presence of nodal lines associated with a d x 2_ y 2 order parameter would be to measure the 
magnetic anisotropy energy, or magnetization torque, for in-plane fields. Consider a velocity field v s — v X ' x' + v y i y'. 
The projections of the velocity along the nodal lines x' and y' imply two occupations of quasiparticle states of differing 
magnitudes, which generate the current, 
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(9) 



The important feature is that the current is not parallel to the velocity field, except for the special directions along 
the nodes or anti-nodes; neither is the magnetic field parallel to the applied field, H. This implies that the magnetic 
energy, U = — g^- J b ■ Hd 3 x, is anisotropic in the ab-plane. For a film in a parallel field, H = H{—sindx' + cosfly'), 
the magnetic anisotropy energy to leading order in H/H a is [ pTTj] ], 

TT2 / TT \ 

U an {-d) = — XA $ {sin 3 tf + cos 3 tf} , < # < tt/2 , (10) 
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where A is the surface area of the film and $ = |(1 + ^cosh(^)) sink 4 (jt)/ cosh 3 (j^), which is constant equal to 
~ 1/3 for d > A, and varies as <3? ~ j(d/2X) A for very thin films. The anisotropy energy is minimized for field directions 
along the nodal lines, and is maximum for fields along the anti-nodes. In order to suppress vortex nucleation, thin 
films with dimensions d < X are desirable; the lower critical field for vortex nucleation is increased by roughly (X/d) 
in a thin film. The optimum geometry might be a superlattice of superconducting/normal layers with an S-layer 
thickness, £ -c d < X. In this case the field at each SN interface is essentially the external field, and the anisotropy 
energy is enhanced by the number of S-layers. 

Torque magnetometry has proven extremely useful for measuring the anisotropy of the penetration lengths for 
current flow along the c-axis compared to the basal plane [ fl5| l. In the linear response limit we do not expect 
anisotropy of the penetration depth for current flow in different directions in the basal plane; however, in-plane 
anisotropy of the current arising from nonlinear field corrections should be observable at low temperatures, if the 
superconducting state has an unconventional order parameter. For comparison, the maximum magnetic torque from 
eq.© is t ~ (1/6V3tt)H 2 (H/H )AX - 10- 3 dyne - cm/rad, for H — H c i = 250 G, A = (2,000 [imf and 
A = 1, 400A, which is small, but comparable to the smallest values of torque reported in [ |15])], 

In summary, we suggest that measurements of the field dependence of the penetration depth for fields parallel to the 
CuO layers at temperatures well below 4 K should provide an answer to the question of whether or not the gap in the 
CuO superconductors has a line of nodes. If so, torque magnetometry and/or the anisotropy of the low temperature 
penetration length could be used to locate the nodal lines on the Fermi surface, thus providing direct evidence for or 
against a d x i_ y i order parameter for the CuO superconductors. 

This research is supported by the Science and Technology Center for Superconductivity through NSF grant number 
DMR 88-09854. 
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Figure Caption 

Figure 1. Phase space for v s \\x' 
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